Horizons Ephemeris Examples

The JPL Horizons platform provides ephemeris data as Cartesian state vectors, osculating orbital elements, and other representations. The precise REST API accepted by the JPL Horizons servers is implemented in Julia with HorizonsAPI.jl.

using HorizonsAPI

fetch_vectors(399; format="text") # Earth's position w.r.t. the solar system barycenter

HTTP.Messages.Response:
"""
HTTP/1.1 200 OK
Server: nginx
Date: Tue, 10 Sep 2024 01:14:39 GMT
Content-Type: text/plain
Transfer-Encoding: chunked
Connection: keep-alive

API VERSION: 1.2
API SOURCE: NASA/JPL Horizons API



*******************************************************************************
Ephemeris / API_USER Mon Sep  9 18:14:39 2024 Pasadena, USA      / Horizons
*******************************************************************************
Target body name: Earth (399)                     {source: DE441}
Center body name: Earth (399)                     {source: DE441}
Center-site name: BODY CENTER
*******************************************************************************
Start time      : A.D. 2022-Jan-01 10:00:00.0000 TDB
Stop  time      : A.D. 2022-Jan-02 00:00:00.0000 TDB
Step-size       : 60 minutes
*******************************************************************************
Center geodetic : 0.0, 0.0, 0.0                   {E-lon(deg),Lat(deg),Alt(km)}
Center cylindric: 0.0, 0.0, 0.0                   {E-lon(deg),Dxy(km),Dz(km)}
Center radii    : 6378.137, 6378.137, 6356.752 km {Equator_a, b, pole_c}
Output units    : KM-S
⋮
9098-byte body
"""

The CENTER parameter specifies the ephemeris reference position.

fetch_vectors(499; format="text", CENTER="@399") # Mars' position w.r.t. Earth

HTTP.Messages.Response:
"""
HTTP/1.1 200 OK
Server: nginx
Date: Tue, 10 Sep 2024 01:14:40 GMT
Content-Type: text/plain
Transfer-Encoding: chunked
Connection: keep-alive

API VERSION: 1.2
API SOURCE: NASA/JPL Horizons API



*******************************************************************************
Ephemeris / API_USER Mon Sep  9 18:14:40 2024 Pasadena, USA      / Horizons
*******************************************************************************
Target body name: Mars (499)                      {source: mar097}
Center body name: Earth (399)                     {source: DE441}
Center-site name: BODY CENTER
*******************************************************************************
Start time      : A.D. 2022-Jan-01 10:00:00.0000 TDB
Stop  time      : A.D. 2022-Jan-02 00:00:00.0000 TDB
Step-size       : 60 minutes
*******************************************************************************
Center geodetic : 0.0, 0.0, 0.0                   {E-lon(deg),Lat(deg),Alt(km)}
Center cylindric: 0.0, 0.0, 0.0                   {E-lon(deg),Dxy(km),Dz(km)}
Center radii    : 6378.137, 6378.137, 6356.752 km {Equator_a, b, pole_c}
Output units    : KM-
⋮
9099-byte body
"""

The response from the Horizons API contains a lot of information! To more simply fetch and inspect ephemeris data from JPL Horizons servers, use the ephemeris method in HorizonsEphemeris.jl.

using Dates, HorizonsEphemeris

ephemeris("earth", now())

(t = [2.460563551862361e6], cal = InlineStrings.String31[" A.D. 2024-Sep-10 01:14:40.9080"], x = [0.9765912607952846], y = [-0.2041721042447016], z = [-0.0883006124611209], ẋ = [0.003445697171843497], ẏ = [0.01534907054386386], ż = [0.006653827452783479])

The output of ephemeris is a NamedTuple, so you can easily pass it along to DataFrames.jl for convenient inspection & data processing. Note the wrt keyword argument as well; it is a simplified interface to the Horizons API CENTER argument.

using DataFrames

data = DataFrame(
    ephemeris("earth", now() - Year(50), now() + Year(5), Day(1); wrt="jupiter")
)